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Abstract - Electrocardiogram (ECG) signal has been 
widely used in cardiac pathology to detect heart disease. The 
ongoing trend of ECG monitoring techniques to become 
more ambulatory and less obtrusive generally comes at the 
expense of decreased signal quality. To enhance this quality, 
consecutive ECG complexes can be averaged triggered on 
the heartbeat, exploiting the quasi-periodicity of the ECG. 
However, this averaging constitutes a trade-off between 
improvement of the SNR and loss of clinically relevant 
physiological signal dynamics. Recent work has attempted 
to utilize wavelet techniques in the analysis of biomedical 
signals including ECG. In this paper Multi- Resolution 
Analysis property of wavelet transform is used to enhance 
ECG signal quality by denoising the signal using Daubechies 
6 Wavelet as best suited mother wavelet at level 4. 

Keywords: Electrocardiogram (ECG), Adaptive Wavelet 
Transform (AWT), Multi- Resolution Analysis (MRA), 
Signal to Noise Ratio (SNR), Daubechies 6 (db 6). 


I. INTRODUCTION 


In ECG, during every heartbeat, P, Q, R, S, T, and U 
waves can be seen. These waveforms result from 
depolarization and repolarization of different parts of the 
heart muscle. The performance of ECG analyzing system 
depends mainly on the accurate and reliable detection of 
the QRS complex, as well as T and P waves. The P-wave 
represents the activation of the upper chambers of the 
heart, the atria, while the QRS complex and T-wave 
represent the excitation of the ventricles or the lower 
chamber of the heart. The detection of the QRS complex is 
the most important task in automatic ECG signal analysis. 
Once the QRS complex has been identified a more 
detailed examination of ECG signal including the heart 
rate, the ST segment etc. can be performed. In the normal 
sinus rhythm (normal state of the heart) the P-R interval is 
in the range of 0.12 to 0.2 seconds. The QRS interval is 
from 0.04 to 0.12 seconds. Normally, the frequency range 
of an ECG signal is of 0.05—100 Hz and its dynamic range 
of 1-10 mV. 


The Q-T interval is less than 0.42 seconds and the 
normal rate of the heart is from 60 to 100 beats per minute. 
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So, from the recorded shape of the ECG, we can say 
whether the heart activity is normal or abnormal. The 
electrocardiogram is a graphic recording or display of the 
time variant voltages produced by the myocardium during 
the cardiac cycle. The P, QRS and T-waves reflect the 
rhythmic electrical depolarization and repolarization of 
the myocardium associated with the contractions of the 
atria and ventricles. This ECG is used clinically in 
diagnosing various abnormalities and conditions 
associated with the heart. The normal value of heart beat 
lies in the range of 60 to 100 beats/minute. A slower rate 
than this is called bradycardia (slow heart rate) and a 
higher rate is called tachycardia (fast heart rate). If the 
cycles are not evenly spaced, an arrhythmia may be 
indicated. If the P-R interval is greater than 0.2 seconds, it 
may suggest blockage of the AV node. The horizontal 
segment of this waveform preceding the P-wave is 
designated as the baseline or the isopotential line. The P- 
wave represents depolarization of the atrial musculature. 
The QRS complex is the combined result of the 
repolarization of the atria and depolarization of the 
ventricles, which occur almost simultaneously. The T- 
wave is the wave of ventricular repolarization, where as 
the U-wave, if present is generally believed to be the result 
of after potentials in the ventricular muscle. So, the 
duration amplitude and morphology of the QRS complex 
is useful in diagnosing cardiac arrhythmias, conduction 
abnormalities, ventricular hypertrophy, myocardial 
infection and other disease states. 


Monitoring and analysis of the ECG has long been 
used in clinical practice. In recent years, the application 
field of ECG monitoring is expanding to areas outside the 
clinic. An example of such an area is at-home monitoring 
of patients with sleep apnea [1]. Also within the clinic, a 
transition in ECG monitoring applications is taking place. 
With developments in sensor technology (e.g., textile 
electrodes and capacitive electrodes), sensors that are 
incorporated in garments or the mattress of an incubator 
{2] have become available. As a result of these new sensor 
technologies, the comfort of the patient is improving 
progressively. Whereas some years ago the patient had to 


reconcile himself or herself with the discomforts of the 
only available technology, nowadays patients prefer the 
more comfortable ways of recording the ECG. However, 
in most cases, this increased comfort comes at the expense 
of signal quality. Electrodes that are incorporated in 
garments generally provide signals with a lower SNR and 
more artifacts than contact electrodes that are glued to the 
body [3]. Another example of ECG signals with a 
typically low SNR is fetal ECG signals, either recorded 
invasively after membrane rupture or noninvasively from 
the maternal abdomen [5]. Some of the SNR and artifact 
problems that arise during these recordings can be 
suppressed by simple, frequency-selective filtering 
[5]-{7]. However, due to the partial overlap of signal and 
noise bandwidths [8], [9], this frequency-selective 
filtering can only help to some extent. Further 
improvement of the ECG can be achieved by exploiting 
its (quasi-)periodicity. Consecutive ECG complexes 
resemble one another and are, moreover, in general 
uncorrelated with noise and artifacts. Hence, by 
averaging several consecutive ECG complexes, the SNR 
can be improved. For additive Gaussian noise, this 
improvement is directly related to the number of ECG 
complexes included in the average [10]. The drawback of 
averaging multiple consecutive ECG complexes is that, 
besides noise, also the physiological dynamics of the 
ECG are suppressed. That is, changes in the ECG that 
originate from physiological events for instance, changes 
in the ST segment that might be associated to metabolic 
acidosis are suppressed in the averaging, complicating 
clinical diagnosis. From this, it is clear that the averaging 
of ECG complexes entails a trade off between the pursued 
increase in SNR and the time scale at which 
physiologically relevant changes in ECG morphology are 
expected to occur. 


Wavelet transform (WT) is a very promising 
technique for time frequency analysis. By decomposing 
signals into elementary building blocks that are well 
localized both in time and frequency, the WT can 
characterize the local regularity of signals [4]. This 
feature can be used to distinguish ECG waves from 
serious noise, artifacts and baseline drift. Using 
aBayesian framework a sequential aVeraging filter is 
developed that, in essence, adaptively varies the number 
of complexes included in the averaging based on the 
characteristics of the ECG signal. The filter has the form 
of an adaptive Kalman filter. The adaptive estimation of 
the process and measurement noise co-variances is 
performed by maximizing the Bayesian evidence 
function of the sequential ECG estimation and by 
exploiting the spatial correlation between several 
simultaneously recorded ECG signals, respectively [13]. 
In this paper, we have used Multi Resolution Analysis 


(MRA) technique of WT to denoise ECG signal and 
enhance signal quality by increasing SNR. 


Il. WAVELET TRANSFORM 
A. Continuous Wavelet Transform (CWT) 


Let x(¢) be a signal defined in L(R) space, which 
denotes a vector space for finite energy signals. R is a real 
continuous number system, Such signals satisfy: 


fls@iarso 


The Wavelet Transform or Continuous Wavelet 
Transform (CWT) of a continuous time domain signal (7) 
is given by X(a,t) which is defined as : 


Forward CWT: 


+0 
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Inverse CWT: 


x(t) = all kan pew (F = dat Q) 
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The analyzing wavelet ‘P(t) term can be written as - 
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Notice that 


Then above equations can be written as: 
Forward CWT: 


X(a,t)= foOv.. *(t)dt (4) 


Inverse CWT: 


xy=* i 


T>-00g=-00 


1 
X(a,T). SW yx (t)da.dt (5) 
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Wherea,t€R and a#Oand 'c'is a constant(0 <¢ 


<0) that depends on the wavelet used. The success of the 
reconstruction depends on this constant called 
admissibility constant, to satisfy the following 
admissibility condition 


c= wo o 
Jy 
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Where ‘P(@) is the Fourier transform and ‘¥*(¢) is the 
complex conjugate of the mother wavelet ‘Y(s), x(¢) the 
signal to be transformed, a and t the dilations (scaling) 
and translations (time-shift) parameters, respectively. 
With a suitable choice of the mother wavelet, the scale 
parameter is proportional to the reciprocal of frequency; 
the translation parameter stands for time [12]. 


A. Discrete Wavelet Transform (DWT) 


The corresponding discrete wavelet transform (DWT) 
of a time domain signal x(f) is given by (j,k) which is 
defined as 


+o = j 
X(jk)= f x(0). v (a ©) 


J j 
co Ja, ay 


The analyzing wavelet ‘Y(t) term can be written as 


1 t—kty ay’ 
Wik (t)= 7 j (7) 
lal % 
Then the expression of Eq. (6) can be written as: 


+00 


X(j,k) = J x(). 5, (tat (8) 
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We can define the DWT of a signal x(t) to be 
the set of analysis coefficients: 


Analysis: C), = J x()y,, (dt (9) 
From these we can recover the signal as: 


Synthesis: x(¢)= >. oan a) (10) 
7 ok 


Assuming existence of a scaling function, @(f) we can 
modify this definition as follows: 

Since the spaces are getting larger’and larger as j goes to 
+oo We can approximate any signal x() closely by 
choosing a large enough value of j = J and projecting the 
signal into V,using the basis 


cA,(m) = J 2(00,,,Odt ay 


co 


From these we can approximately recover the signal as: 


x(t) = Yc4, (mb, (0) (12) 


In effect, we replace the signal x(t) by the approximat 
signal given by the projection coefficients, cA0(m). From 
Eq. (10), (11) and (12), we can write 


xf)= ¥c4,(m).O)» (t) 
= Yc4(h)o, (H+ Ved k Wy () 
k k 


=A,(t)+D, (t) (13) 


As before, we call the signals A(t) and D,(t) the 
approximation and the detail at level-1. We call the 
coefficients cA,(k) and cD,(k) the approximation- 
coefficients and the detail-coefficient at level-1. 


We can further decompose A (t) to get: 


x(t) =A; (t)}* D, (t) 
=Fc4,(k)O, 2. (0+ cD 2(0) 
k k 


+> cD (AW 40 
k 
=A; (t)+ D, (t)+ D, (t) (14) 


We call the signals A,(t) and D,(t) the approximation 
and the detail at level-2. We call the coefficients cA,(k) and 
cD,(k) the approximation-coefficients and the detail- 
coefficients at level-2 [12]. 


C. Multi-resolution Analysis 


Define W, to be set of all signals x(t), which can be 
synthesized from the daughter wavelets y,,(t), 


x)= ¥ x, (15) 


j=—- 


where X,(t)= y CM jx (OD) 
k=-<0 


This leads to various decompositions: 


X(t) = A\(t) + D,(t) 


= 4,(t)+D,(1)+D,(t) 
= A(t) +D,(t)+ D,(t)+ D(t) 
= A(t)+ D,(t)+ D,(t) + D,()+D, (0) 


where D(t), in W,, is called the detail at level i and A(t), is 
called the approximation at level i. 
The approximate coefficients can be computed as 


below: a4. (kK)= (x().0)-14 ) 
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n 


=F c4,(n)(6,,,(0.0,,. 0) (16) 


To complete this calculation we have to compute the inner 
product: 


(6 9(00)-.0(0)) = | ¥2/9(2/t-m2'9(2'1 Kj 


- J W277" (2't—n)o(24't— kat 


(substitute ~ 27"t-k) 
= J V2 0(2s+2k-n)o(s)ds 
(Use the 2-scale equation for O(s) ) 


= i V2 (2s + 2k—n)¥ hy(m)V20(2s —m)ds 


= DA (m) | (25+2k —n)o(2s—m)2ds 
(integral is 0 unless m=n-2k) 

=h,(n-2k) 
The detail coefficients can be computed similarly [12] as 


cD, (k) = (xO j4(D) 


= (Teton, ow ul) 


=¥ cA()(0,, OW) 1,0) 


The Equations (16) and (17) are used for decomposing the 
signal into approximation and detailed coefficients. 


(17) 


III DATA ACQUISITION AND DEMARCATION OF 
INDIVIDUAL ECG COMPLEXES 
To evaluate the used technique AWT, a diversity of 
arrhythmic ECG signals is used. The signals are obtained 


from the Massachusetts Institute of Technology—Beth 
Israel Hospital [11]. To map the AWT's performance as a 
function of the SNR of the recorded signals as well, the 
ECG signals are corrupted with additive Gaussian noise of 
various amplitudes, yielding ECG signals with SNRs 
ranging from —3 to 24 dB. The TWA signals comprise a 
rather ideal dataset for evaluating the performance of the 
AWT. They exhibit relatively high SNR values that can be 
made smaller by additive Gaussian noise, and that 
moreover facilitate quantitative assessment of the 
processed ECG signals. Before defining the individual 
ECG complexes, the QRS complexes need to be detected. 
To facilitate this detection, the SNR of the ECG signals is a 
priori enhanced by linearly combing the signals in such a 
way as to maximize the variance. The linear combination 
with maximum variance is referred to as the principal 
component. The QRS complexes are subsequently 
detected in the principal component as local extrema that 
exceed an adaptive threshold. This adaptive threshold is 
updated continuously by means of an AWT and depends 
on the SNR of the ECG signals complexes in the principal 
component, when the SNR changes, the threshold is 
adapted to prevent noise from exceeding it, in the mean 
time ensuring that the QRS complexes still exceed the 
threshold. 


IV. RESULTS AND DISCUSSIONS 

Noisy signal is shown in figure (1) which is analysed 
with Daubechies 6 Wavelet as best suited mother wavelet 
at level 4 with the help of Wavelet tool of MATLAB 
R2012a. Figure 2 shows approximation signal having 
lower frequency and figure (3) to figure (6) depicts 
detailed signal d,, d,, d,, and d, which contains higher 
frequency component of the signal. From the waveform 
shown in figure 3 it is clear that detailed coefficient dl is 
having higher frequency than that of d,, detailed 
coefficient d, is having higher frequency than that of d, and 
detailed coefficient d, is having higher frequency than that 
of d,. Signal is denoised with Wavelet transform and 
denoised signal is plotted which is shown in figure 7. For 
analysis many mother wavelet like Haar, Coiflet, Symlet, 
Morlet, Mexican Hat, Daubechies 2, Daubechies 4, 
Daubechies 6 and Daubechies 8 have been tried but 
Daubechies 6 (db- 6) provides high SNR. 


Figure | 
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Figure 6 
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V. CONCLUSIONS 


An approach to study ECG signal based on adaptive 
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Figure 7 


ECG signal noise reduction and signal enhancement. 
More tests will be conducted to investigate further its 


wavelet transform is investigated in this paper. Computer 


simulation results show this approach is promising for performance in the future. 
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